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Abstract 

We study two versions of quasicrystal model, both subcases of Jacobi matrices. For 
off-diagonal model, we show an upper bound of dynamical exponent and the norm of the 
transfer matrix. We apply this result to the off-diagonal Fibonacci Hamiltonian and obtain 
a sub-ballistic bound for coupling large enough. In diagonal case, we improve previous lower 
bounds on the fractal box-counting dimension of the spectrum. 



1 Sturmian on and off-diagonal models 

On and off-diagonal models are special cases of the Jacobi operators. Given two real-valued 
sequences a and b, a Jacobi operator H acts on £ 2 (Z) in the following way 

Hip(n) = a(n)ip(n + 1) +a(n - l)ip(n - 1) +6(n)-0(n), n e Z. (1) 

H is associated to a self-adjoint tridiagonal matrix with diagonal entries filled with values of b 
and off-diagonal entries filled with values of a. 

Let (3 G [0, 1] be an irrational number, on and off-diagonal Sturmian models will be denned 
respectively by setting in (OQ) 



a(n) = l, b(n) = \ 1 (l(n + l)0\-[nP\) (2) 

respectively 

a(n) = (Ai - A 2 ) ( [{n + 1)/3J - [nfi\ ) + A 2 , 6(n) = 0, n G Z (3) 

where Ai, A 2 G 1R + and for igK, \_x\ denotes the largest integer smaller than \x\. The term 
Sturmian refers to quasiperiodicity of sequences b in ([2]) and a in ([3]). On-diagonal model is more 
usually named discrete Schrodinger operator and when /3 is the golden mean the Sturmian 

model is oftenly called the Fibonacci Hamiltonian. 

The study of those models was initiated with the introduction of the Fibonacci operator in 
the early 1980's by Kohmoto et al. |KKT] and Ostlund et al. |OPRSS| . At that point in 
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time, the interest in this model was based mainly on the existence of an exact renormalization 
group procedure, the appearance of critical eigenstates and zero measure Cantor spectrum. Only 
shortly thereafter, Shechtman et al. [SBGCJ reported their discovery of structures, now called 
quasicrystals, whose are shown to have aperiodic structure. Sturmian sequences are central models 
of a quasicrystal in one dimension. Indeed, it is aperiodic moreover, it belongs to virtually all 
classes of mathematical models of quasicrystals that have since been proposed. We refer to the 
reader to |BM] for a recent account of the mathematics related to the modelling and study of 
quasicrystals. Thus, the study of Sturmian operators was further motivated by the interest in 
electronic spectra and transport properties of one-dimensional quasicrystals. It is shown that those 
operators give rise to anomalous transport for a large class of irrational number |DTll IM1[ IM2] . 



Moreover, it exhibits a number of interesting phenomena, such as Cantor spectrum of Lebesgue 
measure zero [Slj IS3j IBIST] and purely singular continuous spectral measure |DI4 IDKI4 IBIST) . 
Consequently, apart from the almost Mathieu operator, the Fibonacci operator has been the most 
heavily studied quasi-periodic operator in the last three decades; compare the survey articles 
|D14 ID2[ IS2] . Partly due to the choice of the model in the foundational papers |KKT| OPRSSJ 
the mathematical literature on the Fibonacci operator has so far only considered the diagonal 
model. Given the connection to quasicrystals and hence aperiodic point sets, and particularly 
cut-and-project sets, it is however equally (if not more) natural to study the off-diagonal Sturmian 
model. We do refer the reader to [BMJ for background. For further motivation to study also the 
off-diagonal model, we mention that it has been the object of interest in a number of physics and 
then mathematics recent papers [ELl IEL21 iDahl IDGl IK3T1 [VP]. 

In this paper, we are interested in transport properties and in fractal dimension of the spec- 
trum. More precisely, the reader will find in the following: 

For off-diagonal Sturmian model, 

• We show pseudo spectrum and spectrum are equal (Theorem [1]). 

• We established a link between outside probabilities and transfer matrices in both time- 
averaged (Theorem |2]) and non time-averaged settings (Theorem [3]). 

• In the Fibonacci case, we deduce from previous item a dynamical upper bound for the 
wavepacket spreading (Theorem H]). This bound is sub-ballistic for hopping constants Ai, A2 
well choosen. 

For on-diagonal Sturmian model, 

• We improve previous lower bound for the box-counting dimension of the spectrum QM2J). 
This bound is valid for a large class of irrational number verifying a Lebesgue measure 1 
diophantine condition (Theorem [5]). 

• Considering a smaller class of irrational number, but still Lebesgue measure 1, we are able 
to improve our previous lower bound (Theorem [6]). 
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2 The trace map application 

Denote the continued fraction expansion of (3 by 

= [ax,a 2 , ...] 



a-2- 



and define the best rational approximants, 

p k+ i = a k+1 p k +Pk-i, Po = 1, P-i = 0, 

Qk+i = a k+ iq k + q = 0, g_i = 1. 

A classic tool to investigate one- dimensional model is to write to free equation 

Hip(n) = Ei/j(n) (4) 

where E is a real or a complex number. Equation @ can be rewritten in both on and off-diagonal 
case respectively as 

Vb M = /V(n + 1)\ = ^ - 6(n) -A / ^(n) 

1 J V v>(n) ; v 1 o ; W(n - 1) 



a(n)ip(n) J a(n) \a{n) 2 y ya(n — l)^(n — 1 



Denote T n & (£) = ^ ^ and T£(E) = ^ (? )2 A In all the following, objects 

with exponent a are associated to the off-diagonal model, while those with exponent b are asso- 
ciated to diagonal case. When there is no ambiguity, we will drop it. Then for i — a,b and with 
F^(E) = T^(E) . . .T[(E), one constructs recursively a solution ip % of the equation with 

^(n) = F^(E)¥(0). 

We define the so-called transfer matrices, denoted by M^(E) = F* n (E). One important objet 
is the traces of tranfer matrices, denote x l n = TrM^ and z l n = Tr(M^_ 1 M^). In the following, we 
will drop the index % whenever it is possible for sake of simplicity. 

The transfer matrix A;— evolution follows the simple rule (see e.g. [B IST} URTt |R] ) 

M k = M k _ 2 M* k _ v k > 1. (5) 

/— \ (\ —A 

We can extend this relation to k = by setting Mt x — I ^ Al J and M h _ x - 1 ' 



o f/ 1 v o i 

The evolution of these sequences is derived from taking the trace in (JSJ) and by Cayley-Hamilton 
Theorem (see [BISTJ for details); for k > 

X fc +l = Z k Sa k+1 -l(x k ) - X k -iS ak+1 - 2 (Xk), (6) 
z k+l = Z k S ak+1 {x k ) — X k ^iS ak+1 ^i(x k ), (7) 
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where the Si are Chebyshev second order polynomials 

Si+i(x) = xSi(x) - Si-i(x), VI > 0, S (x) = l, S- 1 (x)=0. 
The initial conditions of these two sequences are given by, (see |DGl \Dah\ IBISTj ) 

a _ A l + A 2 z a ~— (8) 

a£ 1 = 2, 4 = E, 4 = E-\!. (9) 
Sequences x l and z l verify the Fricke-Vogt invariant, namely for k > and E, one has 

(4) 2 + (4 + i) 2 + (4 + i) 2 - 44+i4 + i = (c 4 ) 2 + 4 (io) 

With C a = ( A i-^)(Ai+A 2 ) and c fe = A 
A1A2 

3 Periodic approximants for off-diagonal Sturmian model 

We define the following sequence of sets 

a n , k = {E G K : (TrA^M^ < 2}. 

In case k — 0, those sets are called the periodic approximants of H as cr n o is the spectrum of a 
periodic operator defined as in ([[]) with a g n -periodic sequence a, see \S1\ |BISTJ. 

The pseudo spectrum is defined as 5 ro = {E G K. : {x n (E)} n is bounded } and denote N (E) 
the first index such that \xn (E)\ < 2. 

Lemma 1 7jf|a;o| > 2 and \zq\ > 2, then the two sequences {x n ,z n } n are unbounded. 

Proof. Let us first check the following induction, suppose that at step n, we have \z n \ > 2, 
\x n \ > 2 and \z n \ > |x n _i|. Then, ([7j) implies that for some integer I > 0, 

|^n+l| = \ z nSi(x n ) — X n -iSi-i(x n )\ 

> \ z n\ \Sl{x n ) — Si 

> \z n \{\x n \ -I) 1 . (11) 

Since |x n |, \z n \ > 2, one has \z n+ i\ > \z n \ and \z n+ i\ > \x n \. From (Q, the same argument gives 
|x n+ i| > \z n \ completing the induction. Clearly if the induction is fulfilled, it implies that both 
sequences are increasing and taking logarithme in (ITT]) show they are unbounded. 

We check now the initial condition of the induction. Let first assume that Ai < A2. Then 
|xo = -||| > 2 implies that \zq — > ^ > ^ + ^ = x-%. Now suppose that A2 < Ai, we argue 
that \zi\ > \xq\ and \zi\, \x\\ > 2. First if a\ = 1, then 

Z\ = ZqXq — X-\ 
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and since |xo| > |^-i| and \zq\ > 2, we have that \z\\ > \xo\- Moreover, in that case x\ = zq 
completing the proof. Then suppose that a± = I > 1. The assumption A2 < Ai with |xo|, \zq\ > 2 
implies that \E\ > 2Ai. Then there exists a positive constant e such that E = ±(2Ai + e), 



\zi\ 



\(z x - i-i)5n(io) - z Si- 2 (xo)\ 
E 2 XI 

i S, f /•,.'! — 

Ai 



A1A2 



A1A2 



> 



> 



3Ai_A2^ 
A2 Ai A1A2 
2Ai 

A2 A1A2 



Si- 2 (x ) 



2 Ae 

+ y ) l^-i(^o) - ^-2(^0)1 



2Ai 

^7 



+ c £ (\x \ - 1) 



1-1 



> 2\x \ —2 > \xq\. 

The same argument upon x\ implies \x\\ > 2^ + c £ > 2 completing the proof. 



□ 



Lemma 2 Let x n and z n be the trace map sequence define in 
If there exists an integer N such that 



and 



\,5>0. 



I^at-iI < 2 + 5, \xn\>2 + 5 and \z^\ > 2 + 5, 
then modulus of the two sequences are superexponentially increasing 

\x k +i\ >\z k \> e< 5)Gk ' N + 1 ,Vfc > N, 



:i2) 



with 



Moreover, 



Gk — Gk-i + OfcCfc-2? Gq — 1, G-i — 1. 
Boo = P| 0„,o U (Tn+1,0), any N > N . 



(13) 



n=N 



Proof. First part of the Lemma follows from ©-([7]) (see |M2[ IM1] for a detailed Proof). The 
occurence of 5 plays no role in the result and will be needed for technical reason in section 5. The 
condition upon z^ could also be replaced by |xjv+i| > 2, see [BISTJ for this proof version. 

For the second part, it remains to show is that for E G B^, N (E) is finite, then the result 
follows from the first part of the lemma and identical argument that in the on-diagonal case (see 
|Slt IBIST] ). Suppose that x n is bounded and greater than two in modulus. This impling at any 
rank n, that 2 < |x n +i| < |ar n _i|. Elsewhere, replacing the hypothesis |ar n _i| < 2 < by 
the weaker one |x n _i| < in Proposition 4 in |B~IS T], one can show that the sequence x n is 

growing exponentially. Thus each subsequence on even and odd index has to be decreasing and 
bounded from below by 2, and thus has a limit l\. Then the Fricke-Vogt invariant implies the 
sequence z n also has a limit denoted l 2 - Therefore the point (l 2 , h,h) is a fixed point of the trace 
map and the set of equation coming from T(l 2 , 1%, 1%) = (l 2 , h,h) should be verified. This is easily 
shown to be impossible using the fact that ^(Zi)! > 2 since |/i| > 2. □ 
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Theorem 1 For any X\, A2 > 0, o~(H) the spectrum of H and S ro coincide. 

Proof. Using Lemma [2j one can argue the same as in |S1] , replacing N > by N > Nq anytime 
it appears. □ 

4 Dynamical upper bound for off-diagonal model 

In this section, we are mostly interested in dynamical properties of the wavepacket. It is known 
for this type of aperiodic model that the unitary group ip(t) = e~ ttH ip(0) will spread with time t. 
Here, ip(0) is some well localized initial condition, for example, 5\ the Dirac function on the first 
site of £ 2 (Z). We want to quantify this spreading, thus we recall one usual way to measure it. 

Denote by a(n, t) = \ (e~' ttH 5i, S n )\ 2 the probability for the system to be at the n th site at time 
t. Here {5j}j e z is the canonical basis of £ 2 (Z). 

We denote the outside probabilities by 

P(N,t) = <ri,t), 

\n\>N 

and 

P r (N,t) = <n,t), Pi(N,t) = < n ^) 

n>N n<-N 

Namely, P(N, t) is the probability to be outside the ball of size N at time t. 
For all a G [0, +00], as it is done in [GKT], define 

cw x r • f \nP(t a -2,t) + , lnP(t°-2,t) 

b (a) = — hmmi , b^(a) = — iimsup . 

t^oo \nt t^oo lnt 

The following critical exponents are particular of interest: 

af = sup{a > : S (a) = 0}, = sup{« > : S (a) < 00}. 

They verify < af < a^. In particular, if 7 > a+ then P(t 7 ,t) goes to faster than 
polynomially. af can be interpreted as the (lower and upper) rates of propagation of the essential 
part of ip and a^ as the rates of propagation of the fastest part of ip (see |GKT] ). 

It is also convenient sometimes to consider these definitions in average in time. We define the 
time-averaged probability (a(n,T)) = ^ J °° e~ 2t ^ T a(n,t)dt. Replacing then a(n,t) by (a(n,T)) 
we can then define its time-averaged outside probabilities (P(N, T)) and all the exponents above, 
denoted with a tilda. 

Now main notations are set up, we extent the link between (time-averaged or not) outside 
probabilities and transfer matrices to off-diagonal models and show the following: 

Theorem 2 Suppose H defined in (CQ) with a,b defined in ([3]), and let K > 4 be such that 
o~(H) C [— K+ 1,K — 1]. Then the time-averaged outside probabilities can be bounded from above 
in terms of transfer matrix norms as follows: 

(P{N,T)) <exp{-cN) +T 3 J \\F(n, E + i)|| 2 ^ dE. 



6 



The proof starting point, as in diagonal case (see [DTlJ ), is Parseval formula, 

1 r°° 

(a(n,T)> = — J {({H-E-tf-^SntfdE. 

Denote e = ± and = (H - zl)~ l for z G C\R. Let us assume that T > 1 and therefore 
< £ < 1. Then, we have 

(P r (N,T)) = - / M r (N, E + ie)dE, (14) 

where 

M r (A^) = ^ 1(^)5! A) | 2 = ||x*fl(*)*i|| 2 , 

n>7V 

and Xjv(^) = 0, n < N, XN{n) = I, n > N + I. We need to bound M r (N, z) from above. If the 
energy E is outside the spectrum of H, one can use Combes-Thomas estimate (see e.g. |CT] ). 
Assume that rj =dist(£' + ie, cr(H)) > 1, then we have, 

2 

\(R(z)5i,5 n )\ < - exp(— dmm{dr], l}|n — 1|) 
with some universal positive constant d. Using this estimate, one can see easily that 

M r (N, E + ie)dE < C(K) exp(-dN), d>0, 



L 



E:\E\>K 

and thus the limit goes to for any N(T) = T a , a > 0. The remaining problem is to estimate 
the other part of the integral where E + ie may be e— close to the spectrum of H 

L r (N,e)= I M r (N, E + ie)dE. 

J-K 

We link M r (N, z) to the complex solutions to the stationary equation Hip = zip. 
First consider the truncated operator H^, namely 

Hatpin) = w(n + l)ip(n + 1) + w(n)ip(n — 1) (15) 

where 

w(n) = a(n), n < N, w{n) = A 2 , n > N. 
Denote by Rn{z) the resolvent of the operator Hn and define 

S(N,z) = Hxiv^ivWII 2 . 

Lemma 3 For any E e [-K, K], < e < 1, we have 

M r (N, E + ie)< e- 2 S(N, E + ie), 

where the implicit constant depends only on \i,\2- 
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Proof. Using the revolvent identity, 

R(z)5 1 - R N (z)5 1 = R(z) XN (H - X 2 A)R N (z)5 1 = R(z) X N(dT ±1 )R N (z)5 1 . 

where d is a constant, namely Ai — A2 or and T is the shift operator , Tip(n) = ip(n + 1) acting 
on £ 2 (Z). 

Since < < 1, and the sequence a(n) is bounded, we get 

M r (N,z) < 2S(N,z) + 2\\R(z)( X nH - X 2 A)R N (z)5 1 \\ 2 

< 2S(N, z) + 2e~ 2 \\xn{H - \ 2 A)R N (z)5i\\ 2 

< 2S(N, z) + 2e- 2 \\xN(dT ±1 )R N (z)5 1 \\ 2 

<C(d)e- 2 S{N,z). □ 

The next step is to link the quantity S(N, z) with the solutions to the stationary equation. 
For any complex z, define Uo(n,z) as a solution to Huq = zuo, Uq(0,z) = 0, Uo(l,z) = 1, and 
define Ui(n, z) as a solution to Hu\ = zu%, «i(0, z) = 1, iti(l, z) = 0. Since H and Hn coincide 
for all n < N, their solutions uq, u\ coincide for all n < N + 1. We will consider uq, U\ only for 
such n, thus we will use the same notation for uq, U\ for H and H^. 

Now consider the free equation upon 

X 2 u(n - 1) + X 2 u(n + 1) = zu{n), n>N (16) 

and define 

z±^z 2 - AX 2 

Sl,2 — — 



2A 2 

the roots of the caracteristic polynom of the sequence u(n). Since SiS 2 = 1, one has \si\ > 1 and 
\s 2 \ < 1- Therefore Si, s 2 are the eigenvalues of equation ( 1T6|) whose general solution is 

u(n) = Cis" + c 2 4- 

Lemma 4 Let z = E + is, where E e [— iT, i<T] ; < e < 1. Then, for n > N > 3, we have 

\{R N {z)5 1 ,8 n )\ < 2e- \ ' 21 Jl (17) 

^(^VA^) -n (iV + l,z)| 

and 

|„ ( 7 \\n-N 

Ki^ftA)!^" 1 , Jpfil (18) 

|s 2 (^)wi(A',^) -ui(N + l,z)\ 

Proof. The following formula holds for any off-diagonal model, (see |KKL] ): 

(R(z)5i, 5 n ) = d(z)ua(n, z) + b(z)ui(n, z), n>l, (19) 
(R(z)5i,5 n ) = d(z)u (n, z) + c(z)ui(n, z), n < 1, (20) 

where 

b f g ) = H^l , c{z) = ~ m +^ , d t z ) = -m + (z)m.(z) 

b(m,j r (z) + m-(z)) ' 6(m + (2;) + m_(z)) ' 6(m + (,2) + m^(z)) ' 
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with some complex functions m + (z),m^ (z), called the m-functions which depend on the sequences 
a and b. Since ||i?(z)5i|| < e^ 1 and 

d(z) = (R(z)5 1 ,5 1 ), c(z) = (R(z)5 l ,5 ), 

we get \d(z)\ < e -1 , \c(z)\ < e~ x . Since b(z) = 1 + c(z) and e < 1, we also get \b(z)\ < 2e~ 1 . 
Summarizing, 

\d(z)\<e-\ \c(z)\<e-\ \b(z)\ < 2e~ l . (22) 

One should stress that the bounds (I22p hold for any operator and the constants are universal. 

Consider the operators if at and define = Rn{z)5\ (the vector (f> depends on N, of course). 
Since 

(H N - z)4> = 6i, 

the function 4>(n) = (<f>,8 n ) obeys the equation H^ip{n) = Z(f)(n), n>2. Since w(n) = Xi, n > 
N + 1, <f)(n) obeys the free equation ( TT6l) for n > N + 1. Hence, 

(0(iV + £; + l),0(iV + A;)) T = c 1 s 1 (z) k e 1 + c 2 s 2 (z) k e 2 , k > 0. (23) 

Here e ii2 = (s 1)2 (z), 1) T are two eigenvectors of the matrix corresponding to the equation f[T6"j) . 
and the constants ci, c 2 are defined by 

(0(iV + l),0(iV)f = Cl e 1 + C2 e 2 . 

Since |s 2 (^)| < 1, |si(z)| > 1 and G ^ 2 (Z), the identity ( 1231) implies that Ci = 0, and thus 

(0(iV + 1), 0(iV)) T = Av^M*), if, D N (z) ^ 0. (24) 

On the other hand, (fT9l) implies 

0±(jV + l) = d JV (z)« (iV + l,^+6 J v^)«i(iV + l,^, (25) 

0(iV) = d N (z)u {N,z) + b N {z)u 1 {N,z). (26) 

Remark that d^iz) ^ 0, since it is a Borel transform of the spectral measure corresponding to 
the vector 5\ [KKLJ. The same is true for bx(z), Cn{z) since m + (z),m^(z) are functions with 
positive imaginary part (it also follows from d^iz) ^ and the expressions of 6jv, cat). Using the 
fact that Uo(n + 1, z)ui(n, z) — Uo(n, z)ui(n + 1, z) = 1 for any n (in particular, for n = N), and 
(124")) - (126]) . it is easy to calculate Dn{z): 

D (z) = M^) = MX) (27) 

NK ' s 2 (z)u 1 (N,z)-u 1 (N + l,z) s 2 (z)u (N,z)-u (N + l,zY K ' 



As observed above, the solutions Uo, U\ are the same for H, if n < N + 1. It follows from ( 1231) . 
where c\ = and c 2 = D, that 

(fljv(*)<*iA) = ^Ws 2 ( Z ) n - w , n > iV. 
The result of the lemma follows now directly from (12"2"|) and f l2"Tj) . □ 
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Lemma 5 For z = E + is with E G [-K, K] and < e < 1, /or iV > 3, it/e /lave 

M r (N, z) < C(K)e- 4 ( max ||$(n,z)| 

\3<n<Af 

Proof. The second bound of Lemma |3] and the bound ( TTT1) of Lemma H] yield 

oo 

M r (N, z) < A{K)e- A \s 2 {z)u {N,z)-u {N + l lZ )\- 2 ^\s 2 {z)\ 2k (28) 

fc=0 

< B(K)e- A \s 2 (z)u (N, z) - w (^ + 1, z)\' 2 (29) 

with uniform B(K), since |s2(^)| < 1- Then, one has 

M r (N, z) < C(K)e- A (\u (N, z)\ 2 + \u (N + 1, z)] 2 )' 1 , 

using the fact that Uq is exponentially decreasing in N, the product Uq(N, z)uq(N + l,z) only 
adds some universal constant to C(K). Using f fTB"]) . we can prove a similar bound with u replaced 
by Ui and therefore obtain 

M r (N,z) <C(K)e- 4 \\$(N,z)\\- 2 . 
Since M r (n, z) is decreasing in n, the asserted bound follows. □ 

Proof of Theorem [2], The assertion is an immediate consequence of (fl^|) and Lemma [5] (and 
the analogous results on the left half-line). As sequences a and b in ([3]) are symmetric the same 
can be done for (Pi(N, T)). Note that we can replace [3, N] by [1, N] since this modification only 
changes the i^-dependent constant. □ 

We now show the same result for non time-averaged quantities. As in |DT2j, the key is to 
replace the Parseval formula by a Dunford functional calculus (sometime called Riesz-Dunford) . 

Lemma 6 For every n G Z, t G M, and positively oriented simple closed contour 7 in C that is 
such that the spectum of H lies inside 7, we have 

(e- UH 5 1 ,S n ) = ~^- [ e- Uz ((H-z)-%,S n )dz. 



Proof This is a direct consequence of Dunford functional calculus, see |Dul IDuSl IRSj . □ 

Lemma 7 Suppose H defined in (CQ) and a and b in ([3j) and K > 4 is such that cr(H) C 
[-K + 1,K- 1]. T/ien, 

/•if 

P(iV,t) < exp(-cA^) + 



/ Y^\{{H - E - ir l )- l 8 x ,5 n )\ 2 dE. 

J ~ K n>N 
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Proof. For any t > 0, consider the following contour 7 = 71 U 72 U 73 U 74, where 

ll = {E + iy:Ee [-K, K] , y = r 1 }, l2 = {E + iy:E = -K, ^[-l.f 1 ]}, 

7 3 = {£ + iy:£e y = -1}, ^ = {E + iy:E = K, ye [-l,*" 1 ]}. 

Notice that the spectrum of if lies within this contour and that for 267, Qm(z) < t^ 1 and thus 
| e -«**| < e Lemma [6] then implies, 

l(e-*%A)l<E / K^-z)- 1 *!,^!!^. 

If z G 72U73U74, then we can again apply Combes-Thomas estimates and bound by C exp(-dN). 
The integral over 71 can be estimated using Cauchy-Schwarz inequality: 

^ |<(#-z)-^Mn)IN) <CW j K K \{{H-E-\)- l 8 1 ,5 n )\ 2 dE. □ 

Theorem 3 Suppose H defined in ([T|) and a, b in and K > 4 is such that o~{H) C [-K + 
— 1]. Then the outside probabilities can be bounded from above in terms of transfer matrix 
norms as follows: 

P{N,t) <exp{-dN) +t 4 J (max N \\$(n,E+ |)|| 2 ^ dE. 

Proof. The proof starts by using Lemma [6] and [7] instead of Parseval formula. All the steps are 
then the same that in the proof of the Theorem [2j □ 

5 Application to off-diagonal Fibonacci dynamic 

We apply Theorem 2 and 3 for Fibonacci off-diagonal model and show a non trivial dynamical 
upper bound. Namely we show the following Theorem: 

Theorem 4 Consider the off-diagonal Fibonacci Hamiltonian, that is the operator defined in (CQ) 
and ([3]) with (3 = v ^~ 1 . Assume c = c a > 8, then 

11 — 1 /■ 1 



where £ c = c — 2 + \/ c 2 — 4c + 1 . 

The same holds for non time-averaged dynamical exponent a+. 

Remark 1 Picking c large enough, one obtains a non trivial bound, that is better that the so-called 
ballistic bound 1. 
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Considering Fibonacci special case symplifies greatly the trace map evolution, since one has 
to consider only one sequence of trace as z n = x n+ i. Indeed evolution of the tracemap reduces to 
the following: 

Xn+1 •En-En— 1 %n— 1 

with 

Af + A| E E 

x-i = — — — , x = t-, xx = — . (30) 

Denote by F k the Fibonacci sequence. A direct application of Lemma [1] yields in this case to either 
\%o\ < 2 or < 2 and implies Nq(E) < 1 uniformly in E. This allow to extend description 
in term of periodic band spectrum. We recall some results, classical in on-diagonal case (see 
Appendix B), and extended to off-diagonal Fibonacci case since proof depend mostly on the trace 
map application evolution (see [DahJ). For the same technical reasons that in on-diagonal case, 
we should suppose that c > 4. 

Definition 1 Define a band Bk C <J k fi to be of type A if B k C <Tk-i,o and a band of type B if 
B k C cT/t-2,0- 

This definition exhausts all possibilities as seen in the following lemma: 
Lemma 8 Let c > A, and k > 2. Then 

(i) Each type A band B k C cr^.o contains exactly one type B band B k+2 C o~k+2,o an d 710 other 
bands from erfc+1,0 and o- k +2,o- 

(ii) Each type B band B k C o^ o contains exactly one type A band B k+ i C cr k+lt0 and two type 
B bands from a k+ 2,o positionned around B k+ \ and no other bands from (Tfc+i,o an d CTfc+2,o- 

Proof. The proof inpired from diagonal case |R] can be found in |Dahj . □ 

We recall also some estimates on the band sizes given in [Dah]. Those estimates were there 
used to compute the fractal dimension of the spectrum. 



Lemma 9 Let c > 8, and k > 2. Then, with £ c = c — 2 + y/c 2 — 4c+ 1, we have the following 
inequalities: 

(i) For any type A band B k+ i C 0fc+i,O) E £ -Bfc+i implies 

4+i (E) 



< 2c + 7. 



<{E) 

(ii) For any type B band B k+ 2 C a k+ 2.o,E G B k+ 2 implies 



4 + 2(E) 



< 2(2c + 7). 
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By now, we consider the periodic approximants spectrum in C. 

4 fl = {zeC:\x h (z)\<2 + 6}. 

All the properties keep true replacing 0^0 by a k for some small enough fixed 5 (recall the 
occurence of 5 in Lemma [2]), in particular statement (TTB"]) . A condition on c should be added to 
keep Fricke-Vogt invariant, c> \(5) = [12(1 + 5f + 8(1 + 5f + 4] 1 / 2 . 

The following Proposition states, due to classical Koebe distortion theorem, the height of the 
set o~ k is almost the same that its length. 

Proposition 1 If k > 3, 5 > and c > max(8, X(S)) then there exist constants c$,ds > such 
that 

Fk-l Fk-1 

3=1 j=i 

where {^Ji^xg^! are the zeros of Xf., r k and Rf. are the sizes of respectively the smallest and 
the largest band in a k . 

Proof. The proof is a direct consequence of properties of the functions x k (E) which are proper 
and continuous as polynomials in E and the distorsion theorem of Koebe. See |DT1[ IM2} IM1] for 



details. □ 
We have now all the required tools to finish the proof of the theorem HI 

Proof of Theorem HJ, As x k are real, we have with Proposition [1] 

<y 5 k C {z e C:|3m z\ < R k } C{z6 C:|3m z\ < dF~ l{c) }. 
for a suitable 7(c) . This implies 

4,0 U 4+1,0 ^ G C:|3m z\ < dF^ c) }. (31) 

Let us precise how to choose 7(c). 

From Lemma [8] and Lemma [91 it is easy to bound R k '- 

R k < C fc/2 

We should have R k < dF^ 1 ^ so a suitable 7 can be chosen by taking: 

fclog£ c 



7(c) < limsup 



k 2 log F k 

Remarking that F t behave like leads to 7 (c) < -fe , . 

V J 21og X\± 

For e = 3m z > 0, we get a lower bound for \x n (E + ie)\ uniform in E G [-K, K] C R. For a 
fixed £ > 0, we choose such that dF k 7 ^ V ' < e. With (}3~Tj) . this shows + ie)| > 2 + 5 and 

|xfc+i(i? + ze)| > 2 + 5. As + < 2 + 5 or \x (E + < 2 + 5 from Lemma [U we can 

apply Lemma [2] and thus 

\ Xj \ > e log{1+5)F J- k + 1 Vj > k. 
All this motivates the following definitions: 
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Definition 2 For 5 > 0, T > 1, denote by k(T) the unique integer with 

p7(c) p7(c) 
r k(T)-l < rp < r k(T) 



d d 
and let 

For every v > 0, there is a constant > such that 

N(T) < C u T^viT u . (32) 

Applying Theorem [2] and above estimate, we get 

?k / \ -i 

i3 / / 1 1 »,/,-,.„• s 1 1 2 



P r (iV(T),T)<exp(-ciV(T))+T^ / max : M B (E + ±) | ) dE), 



< exp(-ciV(T)) + T 3 e^ 21og(1+5)F Lv^)J . 



From this bound, it is clear that P r (N(T),T) goes to zero faster than any inverse power of T 
and thus 

5L + 



7(c) 

with z/ arbitrary small. 

One complete the proof without time averaging, using Theorem [3] instead of Theorem [2j □ 



6 Fractal dimension of the on-diagonal model spectrum 

We now consider H defined with (TjQ) and (j2J). Denotes by N(e) the minimal number of balls of 
diameter at most e one need to cover a, then the upper (and lower) box-counting dimension are 
defined respectively by 

diniD cr = hmsup — : , dim R a = hmmt — . 

£ ->o lne lne 

When the limit exists, one denotes it simply by dung (a). 

We give a lower bound of minimal number of balls of some explicit decreasing scale e^, needed 
to cover the spectrum. The first idea to cover the spectrum can be to take into account all the 
bands and take as a scale the smallest band size, but this is a bad idea because this minimal 
length decreases faster than the number of bands grows. A better idea will be to count how many 
bands have a length of order Ek ~ X\ h - This yields to a better lower bound for the box-counting 
dimension of the spectrum associated to a large class of irrational number. 

Theorem 5 Set C h = | £)* =1 log^- +2). We have for any irrational number (3 verifying C = 
limsupCfc < +oo and Ai > 20: 

where a is the spectrum of H. 
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Remark 2 The diophantine condition C < +00 is Lebesgue measure 1 (see Appendix A). In 
particular, it is true for random numbers and degree 2 algebraic numbers. 

As in off-diagonal Fibonacci, one can label bands (see appendix B for precise definitions). The 
following Lemma give precise statement of the counting idea. 

Lemma 10 Denote by n kj i,n kj u and n k ju the number of bands of type I, II and III in 
a kt i, <7fc + i.o, crfe+1,0 with a length greater than e k = 4ro =1 (Ai + 5)~ J/ " ox ~ 4 



For all k, we have the following induction: 

n k +i,i = (ofc+i + l)n ki u + ak+iUkjii, 

Tlk+l,II > l{a fe+1 <2}Wfe,j, 

nk+i,ui = a k +in kt n + (flfc+i — l)nk,m, 

with initial conditions n j = 1, n jj = 0, n jjj = 1. 
Moreover, 

nk,n 7^ Y n k ,in 7^ 0, n k j 7^ 0, n k j > n k ,m, 

and 



n k ,n + n k ju > 



V / 5 + 1 



(34) 

(35) 
(36) 



(37) 



Proof. The induction relation is obvious with Lemma [T21 and Theorem [S] from Appendix B. 
The two first properties are made by induction. Initial conditions give level 0. Assume it is 

true at level n, then as a k +i > 0, n k ji 7^ \j ' n k ju 7^ 0, implies n k+ \j 7^ 0. For the second part, 

if &k+i < 2 then n k+ i t u 7^ 0, else a k+ i > 2 implies n k+ \jn 7^ 0. 

To prove n k j > n k) iii, it suffices to see that n k j > n k j n + n k -i t u + n k -xjn- 

Denote by n k the sum of n k ji and n k jn. For the last property, we argue that for k > 1, 

n k > 2n fc _ 2 + n k - 3 - (38) 

It is easy to verify from (138|) that the sequence n k growth verifies (137|) . 
We now show Using p4|> - (]36|) , we get 

n k ,n = [(a k -i + l)nk-2,n + a k -in k -2,iii] l{a fc <2}> 

n k ,in = (afc - l)(a fc _infe_ 2 ,n + (a fc _i - l)n fe _ 2 ,m) + a k n k „ 2 ,i^{a k -i<2}- 
We distinguish all the cases depending on the values of a k and a k -i in the following table. 



n k > 


Ofc-l = 1 


Ofc-i — 2 


a fc _i = 7 > 2 


a k = 1 


2nA;_.2 + 


3n fc _ 2 + 




a fc = 2 


3n fc _ 2 + 2n fc _ 3 


5n fc _ 2 + 2n fe „ 3 


(27 - ik_ 2 


= A > 2 


(A - l)n k _ 2 + An fc -3 


2(A - l)n k _ 2 + An fc _ 3 


(A-l)(7-l)n fc _ 2 + (A-l)7ijfc_3 



Therefore, the slowest growth for the quantity n k> u + n kj iu is for a k = ak-i = 1 which ends 
the proof. □ 
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Proof of theorem [5] With Lemma [TUl we obtain a bound for Uk,ii + nk,m, that is the number 
of bands of length at least Then by definition of box-counting dimension, we have 

j. +/ wv • t ln ( n k,n + n k , in ) 

dimnicr) > hmmi 

BK ' ~ k -lne fe 

and the stated result. □ 



7 Dimension almost sure 

In this part, we consider only random numbers, that is number with C = 5.04 . . . (see Corrollary 
1 in Appendix A). This class of course is included in the class of the previous section and is still 
Lebesgue measure 1. This first Lemma is purely technical and link the ergodicity of a random 
continued fraction expansion (see Appendix A) with our two step setting. 

Lemma 11 Denote E(X, 7) the event {(02^,^2^+1) = (A, 7)} f or some k > 0. For every natural 

A, 7, we have P(P(A,7)) = - + HX+ ^ 2 )) 2 +7(7+2) - ■ For all e > 0, there exists N big enough, such 
that for every N > N Q , we have the probability: 

P(At least (P(£(A,7)) - e) N couples verify E(\i)) = 1. 

Proof. Denote A(r) the event 

A(r) = {Among the n first there is exactly r couples verifying E(X, 7)}. 



Since d2k and a2k+i are independant random variables, we have P(E(X, 7)) 
and set for convienience y> = P(E(X, f y)), we have 



(ln(2))2 



Then the event {There is less than r couples verifying E(X, 7)} is [Jf, =0 A(k) and 



\fc=0 / k=0 



n—k 



^(^)W£((^TT>) 



k=0 



Using C^ +1 = frfCn — V^T^n compare C% and C r n . Then assuming that r < -|, one shows 
that 

n — r n — n/B n 

< ttt r^4^ r = zz < 1. 



[B- l)(r + l) ~ (S-l)(n/S + l) n + 5 
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It follows that the remaining sum in (139|) is geometric and is bounded by some constant K. 
Using then Stirling formula, one can estimate C r n : 



n" 



[n — r ' r' 



This yields to the following 

P[\jA(k)) <K 



\k=0 



i) \ n f B — 1 \ n /r + 1 



n — r 



B 



If r < (jj — s) n, for some e > 0, then the right term in the above expression tends to 0, which 
ends the proof. □ 

Theorem 6 For almost all irrational number (3, an arbitrary small v > and Ai > 20: 

D-v 



dimj(cr) > 



C + ln(A 1 + 5) 



with D = 1.0382 . 



Remark 3 D is obtained as the limit of an explicit converging series. Value of D is roughly twice 
bigger than the bound (13"3"|) in Theorem^ It is also worthy to be compare it with optimal bound 
for golden and silver mean D = 0.88 . . . |DEGTl ILPWj . 

Proof. First, simplify the recursion relation in Lemma \TU\ by ignoring the smallest term, and 
thus obtain: 

n2k,II + n 2k,III > C2k (n2k-2,II + n2k-2,IIl) 

where C2k depend only on the values of the couple (a2k, Q2fc-i) according to the table: 



C2k{a-2ki tt2k-x) 


a-2k = 1 


d2k — 2 


a-2k = A > 2 


&2k-\ = 1 


2 


3 


A- 1 


&2k-\ = 2 


3 


5 


2(A-1) 


oak-x = 7 > 2 


7 


27-I 


(A-l)( 7 -l) 



Using now Lemma dU and probability value for at recall in Corollary [T] in Appendix A, one 
has for all k > and almost all (3: 



n 2 k 11 + n 2k m > 2 fc ( p ( £ ( 1 ' 1 ))~ E )3' c ( p ( £; ( 2 ' 1 )) +p ( i? ( 1 ' 2 ))- £ )5' c ( p ( i; ( 2 ' 2 ))- £ ) JJ \ k (P(E(hX)-e) 

A=3 

n n n 

J~J(A - i^(P(E(x,i)-e) JJ(2A - i) k ( p ( E ( 2 > x )- £ ) JJ(2A - 2) k(p{E ^' 2) - £) X 



X 



A=3 



A=3 



A=3 



n ((A-i)(/3-i)) fc ( p ^ A )- £ ). 

A,/3=3 
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In the previous equation, we replace all the e constants by only one constant denoted also e 
and we also take the same n in all the products. We can do those simplifications since we consider 
only a finite number of terms. Taking the logarithm on the above expression, one obtains: 

h(5) (ln§)' + ± ( M A) + ln(A - 1)) (h|) In (l + ^L_) + 
± (1„(2A - 1) + h(2A - 2)) (h |) h (l + ^LJ) + 
t(h(A-l) + h W -l))h(l + ^L_)h(l + ^L_ 



\,0=3 

We can rewrite the last double sum: 



(40) 



gg(.n(A - 1) + HP - 1)) in (l + ^) in (l + ^) 
2 gg M A- 1 )h( 1 + ^)h( 1 + ^L_) 



A=3 /3=3 



2£MA-l)h( 1 + ^)gh( 1 + ^L_) 

2 EMA-Dh(i + I ^)(h 



As all the series in n have positive terms and are convergent, we can bound from below by 
taking the limit minus some v > arbitrary small. 

The limit of r.h.s of (HO"]) can be numerically compute and gives D = 1.0382 .... □ 

Appendix A: A theorem by Khintchin 

We recall here some of the ergodic properties of the continued fraction expansion process. 

Theorem 7 |K] Suppose that f(r) is a non-negative function of a natural argumentr, (r=l,2,. . .), 
and suppose that there exist positive contants C and 5 such that 

f(r) < Cr^- & . 

Then, for all numbers in the interval (0, 1), with the exception of a set of measure zero, 

F+2T 



i n 00 In ( 1 + -rr 

^m-Tf{a k ) = Ym K , r(r+ 
»^n^ Jl ; ^ JK ' In 2 

k=l r=l 



18 



Corollary 1 For almost all (3 with respect to Lebesgue measure, 

3 

C = lim sup - log(oj- + 2) = 5.04 . . . 

k K . 

3=1 

Moreover, the density of the number i in the sequence {ak}k is almost surely 

In (l + -r^J) 

Proof. Apply Theorem [7] with f(r) = m(r + 2) and f(r) = lj(r) respectively. □ 

Appendix B: Gap labelling and band estimates 

We recall precise information on periodic on-diagonal spectrum sets. 

Definition 3 For a given k, we call 

type I gap : a band of ct^i included in a band o/er^o and therefore in a gap of o~k+i,o, 
type II band : a band of 0^4-1,0 included in a band of o~k-\ o,nd in a gap of o~k,o, 
type III band : a band of o~k+ifi included in a band ofo~j.fi and in a gap ofo~k,i- 

These definitions exhaust all the possible configurations: 

Lemma 12 [R] At a given level k, 

(i) a type I gap contains an unique type II band of ak+2,o- 

(ii) a type II band contains (c^+i + 1) bands of type I of <Jk+i,i- They are alternated with (afc+i) 
type III bands of ak+2,0 

(iii) a type III band contains (dfc+i) bands of type I of o~j.+i,i- They are alternated with {at+i — 1) 
type III bands of ak+2,0 

Let A = {I, II, III} be a three letters alphabet. For each band B of spectrum at level k, 
correspond an unique word igix . . . ik G A n+1 such that B is a band of type ik included in a band 
of type ik-i at level k — 1,..., included in a band of type io at level 0. This word will be called the 
index of B. More than one band can have the same index. Let T n = {tij(n)) 3 * 3 be a sequence of 
matrix and r = z Zi . . . ik an index, we define: 

L T (T) = ^0^(1)^^3(2) . . ■ti k _ 1 ,i h {k). 

Theorem 8 |LWj If (3 = [0,1,0,2...] is an irrational number in (0, 1) and H defined as in (CQ) and 
(J2} with Ai > 20 then any band B of index t verifies, 

AL T (Q) < \B\ < 4L T (P) 

where P = (P n )„> 

/ c^ _1 \ 

P n = d/On Cx/a n 

\ci/a n cx/a n ) 
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with c\ 



and Q = (Q n ] 



n>0 



Qn 





c 2 (a n + 2) 
,c 2 (a n + 2) 










c 2 (a n + 2)- 3 
c 2 (a n + 2)~ 3 



with c 2 



Ai+5 ' 
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